Properties of atomic intercalated boron nitride K4 type crystals 
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The stability of atomic intercalated boron nitride K4, crystal structures, XBN (X=H, Li, Be, B, 
C, N, O, F, Na, Mg, Al, Si, P, S, CI, K, Ca, Ga, Ge, As, Se, Br, Rb or Sr) is evaluated by the 
geometric optimization and frozen phonon calculations based on the first principles calculations. 
NaBN, MgBN, GaBN, FBN and CIBN are found to be stable. NaBN, GaBN, FBN and CIBN are 
metallic, whereas MgBN is semiconducting. 
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After the study of crystal structure by Wellsi, Sunada 
discovered the 3-D periodic net structure, called a 
crystal, from the graph theory^. Interestingly, the K/^ 
crystal is believed to have three bonds per atom, simi- 
lar to graphite. Therefore, the crystal structure can be 
believed to have features intermediate between those of 
diamond and graphite. However, its existence and prop- 
erties were not investigated in detail. 

Subsequently, Rignanese and Charlier— , and Itoh et 
alA investigated its stability and electronic properties by 
first principles calculations based on the density func- 
tional theorj*^. However, phonon calculations^ii showed 
that this possibility is very small. Although the K/^ type 
crystal structure is thermally unstable, its crystal frame 
is found to be thermally stable when distorted by doping 
impurities such as Na or Mgi. 

Can other type crystal structure be found in nature 
? Mathematically, the 7^4 crystal structure can be con- 
structed by two different types of elements, as shown by 
Kotani et al^. The BN system is the simplest example. 

Crystalline BN is a fascinating material because of its 
high thermal conductivity, extreme hardness, high melt- 
ing temperature and wide band gap. In contrast to dia- 
mond, BN does not react with iron and steel, making it 
a promising candidate for coating high duty toolsets. 

As in the recent situation regarding the carbon and 
atomic intercalated carbon systems, the stability and 
properties K4 crystal for boron nitride and atomic inter- 
calated boron nitride systems have not been investigated. 
Such investigation will provide valuable knowledge about 
the system stability and properties. 

This study examines the structural stability as well as 
several properties of BN and XBN (X = H, Li, Be, B, C, 
N, O, F, Na, Mg, Al, Si, P, S, CI, K, Ca, Ga, Ge, As, Se, 
Br, Rb or Sr) crystals with a A'4 type BN frame by first 
principles calculations. 



First principles calculations based on density func- 
tional theory^ were performed for BN and XBN crystals 
with K4 type of frame using the Vienna Ab-initio Sim- 
ulation Package (VASP)^^. Local density approximation 
(LDA ) M?r ii il?,i is used for the exchange-correlation energy 
functional. In these calculations, all the crystals are con- 
sidered as spin-un polarized systems. 

Computational costs were reduced by employing the 
projector-augmented wave methodic to approximate nu- 
cleus, inner core electrons and valence electrons in each 
atom of the crystal. For H, Li, Be, B, C, N, O, F, Na, Mg, 
Al, Si, P, S, CI, K, Ca, Ga, Ge, As, Se, Br, Rb and Sr, 
valence electrons of ls\ 2s\ 25^, 25^2^1, 25^2^2^ 25^2^3, 
25^2/, 2s22p5, 3s\ 3s2, 8523^1^ 3s^3p^, 3s^3p^, 3s^3p'^, 
Ss^SpS, 3523/43!^ 8523^64^2^ 4s24p\ 4s24p2^ As^Ap^, 
4s24p4^ 4s24p^, 4s24p^5s^ and 4524^65^2 ^^^.^ considered, 
respectively. 

The following procedure is adopted to evaluate the sta- 
bility of the XBN A'4 crystal structure. 

(1) An original primitive unit cell for the BN A'4 crystal 
is shown in Figure 1(a). Four X atoms (X = H, Li, Be, B, 
C, N, O, F, Na, Mg, Al, Si, P, S, CI, K, Ca, Ga, Ge, As, 
Se, Br, Rb or Sr) are allocated to reduced coordinates 

(i h l)' (h i' D' d' h D' (i i D t^e primitive 
unit cell. The primitive unit cell of XBN with P2i3 (T") 
symmetry is shown in Fig. 1(b). The symmetric plane of 
the atomic intercalated system is shown in Fig. 1(c). As 
Fig. 1(d) shows, three adjacent atoms (red) are aligned 
on a straight line at equal distances (B-X and X-N (Ba-X 
and X-Na shown in Fig. 1(e))). 

(2) The binding energy vs volume curve is evaluated 
and fitted by Murnaghan's equation of state^^. For these 
calculations, Brillouin zone integration is performed for 
8x8x8 fc-point meshes generated by the Monkhorst-Pack 
scheme. The method of residual minimization/direct in- 
version in the iterative subspace is used to accelerate 
the convergence of self-consistent total energy calcula- 
tions. The convergence criterion was set to be within 1 
xlO~^ eV/formula unit cell. The cut-off energy for the 
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plane- wave expansion of valence electrons for the primi- 
tive unit cell was determined, so as the number of plane 
waves to be constant over the full range of lattice con- 
stant. Around the minimum of binding energy vs. vol- 
ume curve. The cut off energies were set to be 500 eV 
for BN and XBN crystals. 

(3) The unit cells obtained from step (2) are optimized 
with freedom of the atomic configuration under the sym- 
metry constraint. 

(4) The unit cells obtained from step (3) are optimized 
without any constraint in the crystal structure. The con- 
vergence criterion for step (2-4) was set to be within 1 
xlO"'^ eV/A unit cell. 

(5) Frozen phonon calculations are perfomed using the 
FROPHO codoi^ which is based on the Parlinski-Li- 
Kawazoc methodic. To obtain the force constants for 
the phonon calculations, the atomic displacements are 
set to be 0.01 A . The Born - von Karman boundary 
condition is applied for each obtained primitive unit cell 
with multiplication after step (3) and (4) for the phonon 
calculation. In these calculations, the primitive cells arc 



multiplied by two for each direction of the unit vector. 



III. RESULTS AND DISCUSSIONS 

In the present study, most of the crystal structures 
with a K4 type frame considered for BN and XBN sys- 
tems maintain the type crystal frame composed of B 
and N atoms with distortion through the entire process 
of geometrical optimization. However, SiBN, GeBN and 
RbBN are transformed to stacked BN decagonal ring- 
layer crystal structures with Si, Ge and B-N-B-N-B-N 
strip crystal structures with Rb, respectively, as shown 
in Fig. 2. In all other cases, the crystal symmetry is 
maintained within an accuracy of less than O.OOlA . The 
lattice constant (a), volume at lattice constant (Vb), co- 
hesive energy per atom (Ecoh) and bulk modulus at Vq 
(Bq) are evaluated for the fully optimized structure. In 
this study, Ecoh for XBN, Ecoh.XBN, is defined as 



_ —ExBN + Nxatom X Exatom + NBatom/Ki X EBatom + Nj^atom/Ki X ENatom 

Ecoh.XBN = J- —r —r ■ 

I 



Here, Ea and Na are the total energy and number of 
atoms for A, respectively. To obtain these values, Mur- 
naghan's equation of state is used to fit the binding en- 
ergy vs. volume curve. Parameters in the equation are 
determined using the least squares method in the range 
of about 0.8Vb <V< 1.2Vb, where the root mean square 
is set to be less than 5.0 meV/atom. The obtained values 
of a, Vb, Ecoh and Bq for the crystals are shown in Fig. 
3(a). 

In XBN, the values of a and Vq increase with increas- 
ing atomic number of X. Generally, the values of a for 
XBN crystals are larger than that of boron nitride K4. 
The boron nitride K4 frames are mostly expanded by 
intercalating X to the boron nitride system. 

However, the values of Vb for X = H, Li, .., S, are 
smaller than that of boron nitride K4 whereas those for 
X = CI, .., Sr are larger. This is because boron nitride 
K4 has sufficiently wide vacant spaces that are filled by 
X atoms, and the spaces occupied by X atoms becomes 
larger with increasing atomic number of X. 

On the other hand, Ecoh decreases with increasing pe- 
riod of X. Ecoh is strongly negatively correlated with a 
and Vb, as shown in the correlation coefficients: Ta-E^ah 
= -0.770 and rvo-E,,,, = -0.769. The Ecoh values of XBN 
are smaller than boron nitride K4 values. These results 
suggest that the X-B and X-N bonds in XBN crystals 
are weaker than the B-N bonds in the boron nitride K4 
crystal. For each period of X in XBN, Ecoh generally 



shows the largest values for IIIB(IVB) elements of X and 
the smallest values for I A or VIIB. The maximum energy 
gain for X ~ IIIB(IVB) can be explained by completion 
or over completion of the electronic shell in the X, B and 
N atoms composing the K4 type frame of XBN crystals. 
Subsequently, no correlation appears between Ecoh and 
the dynamical stability. 

Like Ecoh, Bq also decreases with increasing period of 
X. Bq and Ecoh are strongly correlated, as shown in the 
correlation coefficient: rE^^,^-Bo = 0.837. Bq values in 
XBN arc generally smaller than that of BN. 

The trend of each physical value is qualitatively the 
same as in XC2 systems^. However, the correlation coeffi- 
cients for a, Vq, Ecoh and Bq in XBN show slightly smaller 
absolute values than those of XC2 systems: ra~E^ah = " 
0.792, rvo-E,„^ = -0.796 and rs^^^-Bo = 0.848." This 
suggests that the properties of XBN systems can be un- 
derstood by our knowledge of XC2 systems. 

To better understand the correlation between stability 
and structure, the nearest-neighbour distances and an- 
gles between different atoms specified in Figure. 1(e) are 
investigated [d B-N, d X-Ba, d X-Bs, d X-Bc, d X-B 
(the average of all d X-B values), d X-N^, d X-Ns, d 
X-Nc, d X-N (the average of all d X-N values) and d X- 
X, and the angles Z X-B^-Nc, Z X-N^i-Bc and dihedral 
angle Z B-N-B-N]. 

Surprisingly, the relative position of X in the crystal is 
maintained, yielding Va-dx-x = 0.999. That is, d X-X 
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determines the volume of the crystaL Therefore, d X-X 
can be considered as a standard distance in XBN systems. 
To emphasise the similarity of the structures, the various 
distance ratios of the nearest-neighbour distances to the 
d X-X are shown in Figs. 4(a) and (b). 

As Figure 3(b) shows, d B-N, d X-B, d X-N and d X-X, 
like a and Vq , show certain trends with increasing atomic 
number of X. d B-N, d X-B and d X-N are very strongly 
correlated with a, as shown in the correlation coefficients: 
Ta-dB^N = 0.842, r^_ax-i3 = 0.973 and r^_2x-N = 
0.990. Based on the values of d X-X, the values of d X-B 
and d X-N are possibly determined, d B-N seems to be 
determined indirectly by d X-B, d X-N and d X-X. On 
the other hand, these distances have weaker correlations 
with Ecoh- rE,,f,-dB-N = -0.592, r^^^^_2x^B = -0-699, 
and ^-'dX-N ~ -0.723. Also note that d X-B and 
d X-N, like a and unlike d B-N, show strong correlation 
with Ecoh- 

Here, the unequal atomic distances between X and 
B atoms and X and N atoms (d X-B^, d X-B^, 
d X-Bc, d X-Na, d X-Ns and d X-Nc) should be 
noted. It is apparent that the correlations of d X- 
B and d X-N with a and Ecoh result from the cor- 
relation of d X-Ba, d X-Bb, d X-Bc, d X-N^, d X- 
Ns and d X-Nc with a and Ecoh as shown in the 
correlation coefficients: ra-dX-BA — 0.543, Va-dX-Bs 
= -0.042, ra-dX-Bc = 0.927, Va-dX-N^ = 0.812, 

ra-dX-Ns = 0.338, Ta-dX-Nc ~ 0.864, rE^„,,-dX-BA 

= -0.436, rE^^^-dX-Bg = 0.161, rE^„^-dX-Bc = -0.781, 
rE,„h-dX-NA = -0.808, TE^^^-dX-NB = -0.073, and 
rE^^^-dX-Nc = -0.629. These values indicate that d X- 
N^, d X-Bc and d X-Nc {d X-B^) are mainly determined 
by the volume and the energy gain. 

For each group of X, d X-B a, d X-B^, d X-Bc, d X- 
N^, d X-Ns and d X-Nc increase with increasing atomic 
number of X. Here, for each period of X, d X-B^ (d X-N^) 
increases(decreases) up to X = IVB or VB and then de- 
creases (increases). Apparently, this trend is determined 
by the number of valence electrons of X. However, an 
opposite trend is shown in a system with X of second 
period. This can be attributed to the similarity between 
elements of K4 type frame and X. Here, we should note 
the strong (weak) negative correlation between d X-B^i 
and d X-Bb {d X-N^ and d X-Ns): rdx-BA-dX-Bs = 
-0.818 {rdx-NA-dX-Ns 

= -0.134). 

As Table HI in the Appendix shows, for each third 
and forth period with increasing atomic number of X, 
the dihedral angle Z B-N-B-N decreases up to X = 
IVB and then increases. For each period, the dihe- 
dral angle becomes minimum at X = IVB, which breaks 
the Ki type frame structure. For the values of d X- 
Ba, d X-Bb, d X-N^ and d X-N^ which have larger 
amplitudes compared with d X-Bc and d X-Nc, the 
correlation with Z B-N-B-N are investigated. Direct 
correlation between d X-B^, d X-B^, d X-N^, d X- 
Ns and Z B-N-B-N: rax-BA-^B-N-B-N = -0.798, 

TdX-Bs-^B-N-B-N — —0.821, TdX-NA-^B-N-B-N — 



0.196, and rdx-Ns-^B-N-B-N = —0.145 shows much 
preference of d X-B^ and d X-B^ than d X-N^ and 
d X-Nb. X dependence of the angle Z X-B^-Nc and 
Z X-N^-Bc have stronger correlation with dihedral an- 
gle Z B-N-B-N: t^x-Ba-Nc-^b-n-b-n = 0.888 and 
r^x-NA-Bc-^B-N-B-N = -0.662. Z X-Ba-Nc(Z X- 
Nyi-Bc) depends strongly on that of d X-B^ and d X-N^i: 

TdX-BA-^x-BA-Nc = -0.964 and rdx-NA-zx-NA-Bc 
= -0.772. These suggest that X dependence of d X-N^ 
is indirectly correlated to Z B-N-B-N, through and Z X- 
N^-Bc while d X-Ba is directly correlated. The values 
are determined by the number of outer-most valence elec- 
trons and the occupied space of X. Weak correlation of Z 
B-N-B-N with a, Vq, Ecoh and Bq should be also noted. 

The distinct difference between X = lA and HA and 
the other X atoms found in XC2 systems with X of the 
third, forth(fifth) periods in the periodic tabled are not 
found in XBN systems. Note the exceptional features of 
BeBN in the distance and angle trends. Apparently, this 
results from the higher complexity of the bonds formed in 
XBN systems relative to those of XC2 systems. However, 
the relative atomic configurations of XBN systems can be 
considered to have some deviation from XC2 on the basis 
of a comparison between the normalized distances shown 
in Fig. 2(c) in ref<i and Figs. 4(a) and (b). 

In Figure 5, the phonon density of states(DOS) with- 
out step (4) for BN and XBN (X = H, Li, Be, B, C, N, 
O, F, Na, Mg, Al, Si, P, S, CI, K, Ca, Ga, Ge, As, Se, Br, 
Rb or Sr), with K4 type frame of B and N, are depicted. 
In general, the bandwidth of XBN depends on the pe- 
riod of X, which is the narrowest in third period. Thus, 
the inter-atomic distances approximately determine the 
vibrational frequency range in these crystals. As shown 
in the figure, imaginary frequency modes rarely appear 
in NaBN, MgBN, GaBN, FBN and CIBN. 

For the obtained structures in BN, NaBN, MgBN, 
GaBN, FBN and CIBN after step (4), the DOS and dis- 
persion relationship are also evaluated. As Fig. 6 shows, 
the phonon DOS for the fully optimized structures of 
other XBN crystals show almost the same shapes as those 
of the optimized structures from step (3). Therefore, as 
in the previous studyi, the discussion in the previous 
paragraph is expected to be significant, and omitting step 
(4) is effective for reducing the computational cost. How- 
ever, exceptional cases such as SiBN, GeBN and RbBN 
exist, as shown in Fig. 2. 

Figure 6 shows that boron nitride K4 has imaginary 
frequency modes at a wide range of k points. This sug- 
gests a short life-span of the boron nitride K4 in nature, 
similar to carbon K^^. On the other hand, although 
NaBN, MgBN, GaBN, FBN and CIBN have imaginary 
frequencies in the acoustic modes around the F point, no 
imaginary modes appear at the other k points. This sug- 
gests that the structures are stable, although they break 
against phonon vibration with a long wavelength limit. 
From previous calculations for XC2 and SrSi2^, the stable 
existence of NaBN, MgBN, GaBN, FBN and CIBN K4, 
type crystals in nature is expected although maximum 
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values of Ecoh are not realized for the considered XBN 
crystals. 

From the maximum frequencies obtained for BN di- 
amond (1300 cm^^ii) and graphite (1600 cm"^^), the 
orders of the B-N chemical bonds in XBN are consid- 
ered as at most 1. Therefore, the orders are expected 
to be less than 1 for NaBN, MgBN, GaBN, FBN and 
CIBN crystals. The feature of the chemical bonds and 
the stability of K4 systems are analyzed by examining 
the valence charge distribution. 

Figure 7 shows the various type of distribution in 
the primitive unit cells for fully optimized BN, NaBN, 
MgBN, GaBN, FBN and CIBN crystals with a distorted 
K4 type BN frame. Fig. 7(a) details the planes selected 
to show the distributions of the quantities noted in Figs. 
7(c), (e) and (g). In Figs. 7(c) and (e), the high and low 
populations are neglected for better comparison. 

Figures 7(b) and (c) show the isosurface and contour 
of charge density. Apparently the valence charge is lo- 
cally distributed along the lines between adjacent B and 
N atoms for the considered systems. The charge is lo- 
calized more to the N atoms than B atoms because of 
the electronegativity difference. These features are also 
indicated in the zinc blende BN crystalsi^. Although sig- 
nificant charge accumulation is shown for X in FBN and 
CIBN, the difference in the charge distribution between 
BN and XBN is not so large in the BN K4 frame. 

Figures 7(d) and (e) show the differences between the 
valence charge density of the crystals (ex. PNaBN{i^)) and 
their separated components: the induced impurities (ex. 
/97Va(r)) and the frame (ex. pB7v(r))- Let us examine 
the distribution around the X atoms. Significant excess 
charge accumulates at the X atom only in FBN and CIBN 
and not in NaBN, MgBN and GaBN. Between adjacent 
X atoms the charge excess is low. The higher excess and 
depressed charge for FBN compared to CIBN seem to be 
related to the wider phonon frequency range. From this 
example, it can be expected that a higher ionicity of X 
makes the system more robust. 

Charge excess regions appear on the line between ad- 
jacent X and atoms, and the region is closer to the 
X atoms for systems with larger electronegativity of X: 
FBN, CIBN, GaBN, MgBN and NaBN. Further, as shown 
in these figures, significant excess charge appears above 
the opposite side of the B^ atom in the pseudoplane de- 
fined by the B atom and three adjacent nearest-neighbour 
N atoms in NaBN, MgBN, GaBN, FBN and CIBN. The 
resulting significant charge transfer causes ionic bonds 
between X and B^^. For these cases, a region of depressed 
charge exists between the X and BN K4 type frame. The 
charge differences between the X and N^ atoms do not 
show much excess for NaBN, MgBN and GaBN, and are 
depressed for FBN and CIBN. 

The charge accumulations in the region of the tr-type 
orbitals between adjacent B and N atoms are depressed 
for NaBN, MgBN and GaBN but in excess around the 
region close to the N atom for FBN and CIBN. However, 
features of the 7r-type bonds between adjacent B and 



N atoms are strengthened in NaBN, MgBN and GaBN, 
whereas these bonds are not found in FBN and CIBN. 
A correlation exists between the values of the dihedral 
angle Z B-N-B-N and the tt bonds. 

The nature of the bonds is also investigated using the 
electronic localization function(ELF). Figs. 7(f) and (g) 
show the isosurface (0.8 [-]) and the contour. The ELF 
distributions of BN, NaBN, MgBN, GaBN, FBN and 
CIBN have some similarities despite of being significantly 
different. The ELF accumulates distinctively around for 
Ga, F and CI atoms but not so much for the Na and 
Mg atoms in XBNs. First, let us note the distribution 
between the X atom and adjacent B^i and N^ atoms. A 
region of significant ELF concentration appears between 
adjacent X and B^ atoms for NaBN, MgBN, GaBN, FBN 
and CIBN. Between adjacent X and N^ atoms, a re- 
gion of concentrated ELF appears more closer to N^ in 
NaBN, MgBN and GaBN and more closer to X in FBN 
and CIBN. Between adjacent B and N atoms, features of 
the weakened a and strengthened tt character appear for 
NaBN, MgBN and GaBN as expected in XC2 systemsl, 
whereas only the a character appears for FBN and CIBN. 

Valence electron donation from X to B^ in NaBN, 
MgBN and GaBN and from B^ to X in FBN and CIBN 
result in different types of ionic bonds between them. 
The difference seems to originate mainly from the num- 
ber of outermost valence electrons and the order of elec- 
tronegativity in X. d X-B^ becomes larger for electron 
donation from X in NaBN, MgBN and GaBN. However, 
d X-Ba of FBN and CIBN is smaller than the other com- 
pounds. We should also note that d X-B^/d X-X of FBN 
and CIBN is smaller than that of BN, as shown in Figs. 
4(a) and (b). The X = F and CI forms an electronic 
shell, which results to strengthen tt bonds and weaken a 
bonds between adjacent B and N atoms through reduc- 
ing d X-B^ and increace Z B-N-B-N. As shown here, the 
bonds between X and N^ atoms are different from those 
between X and B^i atoms. The greater number of elec- 
trons in Ga compared with NaBN and MgBN shortens 
the relative distance between the X and N^ atoms. This 
results in the formation of the smallest Z B-N-B-N and 
the distinctive tt bonds feature between adjacent B and 
N atoms among NaBN, MgBN, GaBN, FBN and CIBN. 
The nature of the bonds in those crystals are gradually 
different with the number of valence electrons in X. Here, 
we should note again that the net bond orders of XBN 
arc less than that of BN from the phonon calculations. 

Let us discuss the reason for the instability and stabil- 
ity of the BN and XBN systems. In the BN crystal, 
the incomplete tt bonds make the structure thermally un- 
stable as in the C K4 crystal. The thermal stability of 
NaBN and MgBN can be attributed to the same reason as 
that of NaC2 and MgC^i. Valence electrons are donated 
from Na or Mg to the BN K4 type frame. In the K4 type 
frame, polarized bonds between adjacent B and N atoms 
are changed to strengthen the tt bonds and weaken the a- 
type bonds. The crystal structures become stable for the 
formation of these tt bonds. Unlike NaBN and MgBN, 
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KBN and CaBN are thermally unstable. The difference 
can also be attributed to the amount of space occupied 
by the dopant X in XBN. 

The thermal stability of GaBN can be understood by 
comparing it with other systems. The instability of the 
subsequent crystal AlBN is believed to originate from 
the stronger electronegativity of Al compared to Na and 
Mg. In AlBN, Z B-N-B-N is very smaU for the strong 
repulsion of B^ by adjacent X. The loss of the bond- 
ing balance puts the structure on the saddle with highly 
degenerated negative curvature on the potential surface. 
Ga is less electronegative than Al, which results in less 
B^ repulsion. The larger occupied space of Ga in GaBN 
compared with that of Al in AlBN is attributed to the 
greater amount of electrons in the inner core of Ga than 
in Al. These differences seem to retain the K4 type frame 
structure in GaBN and make the Z B-N-B-N value the 
same as that of MgBN. The system seems to be stabilized 
by not only the 7r-type bonds between adjacent B and N 
atoms but also the Ga-N^i bonds. Note that the large d 
X-Ba value of GaBN compared with d X-N^ which is as 
same as that of d X-Ca of Ga.C2^. 

In FBN and CIBN, the F and CI distinctively attract 
B^ and repel Na atoms. This is due to the strong elec- 
tronegativity of X, which attracts the valence electrons of 
neighbour B^, this results in the closing of the X atom's 
electronic shell. The F atom attracts electrons more than 
the CI atom because of its larger electronegativity, and 
this results in smaller d X-Ba and larger Z B-N-B-N 
in FBN. Between adjacent B and N atoms, although tt 
bonds do not form, the a bonds are strengthened. In 
FBN and CIBN, the X-B^ bonds and a bonds between 
adjacent B and N atoms probably stabilize the crystal 
structures. 

The stabilization of the atomic intercalated BN K4 
crystals can be explained by the delicate balance of the 
bonding. Discriminative stability seems to be realized 
by the proper atomic radii of X, although the interpreta- 
tion of thermal stability of the crystal structure is more 
difficult than in the case of XC2 systems. Future accu- 
mulation of knowledge will elucidate the stability of these 
systems. 

Figure 8 shows X-ray diffraction (XRD) patterns of 
the fully optimized crystal structure of the BN, NaBN, 
MgBN, GaBN, FBN and CIBN. Monochromatic radia- 
tion with wave length 1.541 A is assumed in these cal- 
culations. As shown in this figure, the compounds have 
specific XRD peaks at larger 29 than pure BN crystals. 

Figure 9 shows the electron band structures, DOS, and 
local DOS of valence electrons for the atoms composing 
the K4 type frame and the intercalated atoms in the fully 
optimized BN, NaBN, MgBN, GaBN, FBN and CIBN 
crystals. From the bottom energy levels to the upper 
levels, the angular momentum for the electronic states 
changes successively in the order of s, p and d character, 
although every component appears in most of the states. 
From the bottom energy levels to those around the Fermi 
level, the band structures are similar in these crystals. 



NaBN, GaBN, FBN and CIBN K4 crystals are metallic, 
whereas MgBN is semiconducting. 

IV. CONCLUSIONS AND SUMMARY 

The stability of atomic intercalated boron nitride K4 
type crystals, XBN (X==H, Li, Be, B, C, N, O, F, Na, 
Mg, Al, Si, P, S, CI, K, Ga, Ga, Ge, As, Sc, Br, Rb 
and Sr), are evaluated using geometric optimization and 
frozen phonon analysis based on first principles calcula- 
tions. The phonon calculations show that NaBN, MgBN, 
GaBN, FBN and CIBN are stable. NaBN, GaBN, FBN 
and CIBN are metallic while MgBN is semi-conducting. 

V. APPENDIX 

The values presented in Figures 3 and 4 are listed in 
Table I, II, and III. Distance between adjacent B^, X and 
Na atoms for fully optimized BN, NaBN, MgBN, GaBN, 
FBN and CIBN are shown in Fig. 10. 
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FIG. 1: (a) Primitive unit cell of the boron nitride (BN) K4, crystal with P4i32(0^) symmetry, (b) Initial structure of the 
primitive unit cell of the XBN crystal with Ka frame of BN with P2i3(T*) symmetry, (c) The 2x2x2 super-cell of (b) showing 
symmetric plane of the impurity intercalated system, (d) The three adjacent atoms (B-X-N) aligned in a straight line at equal 
distances (B-X and X-N) in (c) are encircled with a thick red line, (e) One of the extracted structure. The B and N atoms can 
be classified into Ba, Bs and Be, and Na, N_b and Nc. 
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SiBN GeBN RbBN 



FIG. 2: 2 X 2 X 2 super-cells of fully relaxed SiBN, GeBN and RbBN from K4, type crystal structures. 



TABLE 1: Determined a (lattice constant), Vq (volume at a), Ecoh (cohesive energy) and Bq (bulk modulus at Vo) for the fully 
optimized crystal structures within LDA. 



Component 


a [A] 


Vo [A Vatom] 


Ecoh [eV/atom] 


Bo [10- ^] 




4.082 


8.502 


7.702 


2.67 


BN 


4.121 


8.762 


7.370 


2.35 


HBN 


4.106 


5.767 


5.532 


2.52 


LiBN 


4.330 


6.766 


5.539 


2.08 


BeBN 


4.287 


6.565 


6.036 


2.50 


BBN 


4.388 


7.042 


6.736 


2.61 


CBN 


4.366 


6.936 


6.591 


2.87 


NBN 


4.372 


6.962 


6.139 


3.02 


OBN 


4.281 


6.539 


6.468 


3.15 


FBN 


4.384 


7.053 


6.144 


2.53 


NaBN 


4.596 


8.109 


4.843 


1.78 


MgBN 


4.631 


8.291 


5.200 


1.86 


AIBN 


4.671 


8.491 


5.645 


1.97 


SiBN (Not Ki type) 




9.515 


6.452 


2.06 


PBN 


4.687 


8.580 


5.447 


2.29 


SBN 


4.642 


8.336 


5.394 


2.34 


CIBN 


4.794 


9.202 


4.798 


1.91 


KBN 


4.986 


10.327 


4.114 


1.17 


CaBN 


4.924 


9.946 


5.317 


1.40 


GaBN 


4.722 


9.087 


5.326 


2.04 


GeBN (Not Ka type) 




10.916 


6.213 


1.73 


AsBN 


4.840 


9.448 


5.041 


1.99 


SeBN 


4.891 


9.753 


4.559 


1.90 


BrBN 


4.987 


10.336 


4.346 


1.71 


RbBN (Not A'4 type) 




20.057 


5.103 


0.65 


SrBN 


5.220 


11.855 


4.844 


1.10 
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FIG. 3: (a) Determined lattice constant(a), volume at a{Vo), cohesive energy(_Ecoh), and bulk modulus at Vo{Bo) (b) and (c) 
Nearest-neighbour distances and angles for the fully optimized B, BN and XBN crystal structures with a K4 type frame within 
LDA. Classification of the B and N atoms (Ba, Bs and Be, and Na, Ns and Nc) is based on the irreducibility of the crystal 
symmetry as shown in Fig. 1. 
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FIG. 4: (a) and (b) Distance ratios of the nearest-neighbour distances for fully optimized BN and XBN crystal structures with 
K4, type frame within LDA. Classification of the B and N atoms (Ba, Bs and Be, and Na, Ns and No) is based on the 
irreducibility of the crystal symmetry as shown in Fig. 1. 
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FIG. 5: Phonon density of states(DOS) of XBN (X=H, Li, Be, B, C, N, O, F, Na, Mg, Al, Si, P, S, CI, K, Ca, Ga, Ge, As, Se, 
Br, Rb and Sr), BN and C crystals with K4, type frame of BN and C. The structures are obtained from geometric optimization 
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FIG. 6: (Left) Phonon DOS of partially and fully optimized structures of BN, NaBN, MgBN, GaBN, FBN and CIBN with a 
Ki type frame of BN. (Right) Phonon dispersion relationship for the fully optimized structures. 
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FIG. 7: (a) Primitive unit cells for fully optimized BN, NaBN, MgBN, GaBN, FBN and CIBN crystals with a K4 type BN 
frame. The planes are selected to show the distribution of the valence charge density, differences in the valence charge density of 
the crystals (ex. pjvaS]v(r)) and their separated components: the induced impurities (ex. pNa{r)) and the frame (ex. pfl]v(r)), 
and the electronic localization function (ELF). Isosurfaces (b) and contours (c) of the valence charge density (1.0 e/A^) for 
BN, NaBN, MgBN, GaBN, FBN and CIBN. The isosurfaces (d) and contours (e) of the differences in valence charge density 
(+(-)0.05e/A^ coloured red(blue)). Isosurfaces (f) and contours (g) of the ELF (0.8 [-]). In (c) and (e), the high and low 
populations are neglected for comparison at the same density scale. XCrySDeiji^ was used for the visualization. 
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FIG. 8: XRD patterns for the fully optimized structures of BN, NaBN, MgBN, GaBN, FBN and CIBN with K4 type BN 
frames. 
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FIG. 9: Electron density of states and the dispersion relationsliip for the fully optimized primitive unit cells of BN, NaBN, 
MgBN, GaBN, FBN and CIBN with type BN frames. 
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TABLE II: Determined nearest-neighbour distances and angles for the fully optimized C, BN and XBN crystal structures with 
type frame within LDA. Classification of the B and N (C) atoms (Ba, Bb and Be, and N^i, Ns and Nc (Ca, Cs and 
Cc)) is based on the irreducibility of the crystal symmetry, as shown in Fig. 1. For BN (C), X means the atomic position 
intermediate to the nearest Ba and Na (Ca) atoms. 



Component 


d B-N [A ] 


d X-B [A ] 


d X-Ba [A ] 


d X-Bs [A ] 


d X-Bc [A ] 


d X-X [A ] 


C 


1.443 


2.209 


1.768 


2.282 


2.282 


2.500 


BN 


1.457 


2.230 


1.784 


2.304 


2.304 


2.523 


HBN 


1.516 


2.267 


1.381 


2.574 


2.256 


2.514 


LiBN 


1.553 


2.318 


2.282 


2.341 


2.308 


2.675 


BeBN 


1.559 


2.295 


2.551 


2.091 


2.413 


2.645 


BBN 


1.617 


2.431 


1.635 


2.909 


2.218 


2.708 


CBN 


1.637 


2.431 


1.497 


2.918 


2.256 


2.682 


NBN 


1.678 


2.455 


1.357 


2.985 


2.292 


2.679 


OBN 


1.612 


2.391 


1.373 


2.863 


2.258 


2.622 


FBN 


1.584 


2.423 


1.396 


2.640 


2.548 


2.705 


NaBN 


1.642 


2.469 


2.269 


2.461 


2.543 


2.820 


MgBN 


1.687 


2.486 


2.498 


2.289 


2.680 


2.838 


AlBN 


1.730 


2.500 


2.791 


2.280 


2.622 


2.883 


SiBN (Not Ki type) 


PBN 


1.906 


2.535 


3.088 


2.064 


2.821 


2.885 


SBN 


1.806 


2.496 


2.799 


2.150 


2.742 


2.849 


CIBN 


1.706 


2.617 


1.811 


2.771 


2.732 


2.941 


KBN 


1.778 


2.682 


2.381 


2.634 


2.831 


3.051 


CaBN 


1.779 


2.643 


2.552 


2.471 


2.846 


3.013 


GaBN 


1.712 


2.537 


2.474 


2.367 


2.728 


2.893 


GeBN (Not Ki type) 


AsBN 


1.789 


2.601 


2.571 


2.373 


2.839 


2.962 


SeBN 


1.742 


2.657 


1.984 


2.761 


2.777 


2.997 


BrBN 


1.769 


2.712 


1.968 


2.833 


2.840 


3.056 


RbBN (Not A'4 type) 


SrBN 


1.885 


2.804 


2.699 


2.629 


3.013 


3.196 
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TABLE III: Determined nearest-neighbour distances and angles for the fully optimized C, BN and XBN crystal structures 
with K4 type frame within LDA. Classification of the B and N atoms (B^, Bs and Bo, and N^, Ng and Nc) is based on the 
irreducibility of the crystal symmetry as shown in Fig. 1. 



Component 


d X-N [A ] 


d X-N A [A ] 


d X-Ns [A ] 


d X-Nc [A ] 


Z X-B.4-Nc n 


Z X-Na-Bc [°] 


Z B-N-B-N [°] 


C^ 


2.209 


1.768 


2.282 


2.282 


90.000 


90.000 


70.529 


BN 


2.230 


1.784 


2.304 


2.304 


90.000 


90.000 


70.529 


HBN 


2.270 


1.415 


2.625 


2.199 


98.641 


100.559 


71.233 


LiBN 


2.355 


1.909 


2.186 


2.673 


86.183 


82.908 


75.318 


BeBN 


2.342 


1.612 


2.385 


2.543 


71.439 


97.469 


30.344 


BBN 


2.418 


1.437 


2.594 


2.569 


104.335 


92.939 


82.667 


CBN 


2.419 


1.358 


2.713 


2.478 


104.418 


97.307 


76.571 


NBN 


2.435 


1.306 


2.813 


2.434 


106.143 


99.607 


71.797 


OBN 


2.366 


1.394 


2.679 


2.378 


105.311 


97.084 


75.843 


FBN 


2.355 


2.174 


2.455 


2.315 


101.720 


83.780 


99.889 


NaBN 


2.485 


2.108 


2.381 


2.715 


86.340 


84.430 


72.238 


MgBN 


2.505 


2.021 


2.511 


2.661 


76.250 


92.110 


44.435 


AlBN 


2.545 


1.876 


2.497 


2.815 


72.778 


93.194 


37.934 


SiBN (Not K4 type) 


PBN 


2.568 


1.670 


2.709 


2.726 


60.840 


103.998 


12.373 


SBN 


2.541 


1.684 


2.722 


2.645 


66.087 


103.532 


17.301 


CIBN 


2.579 


2.276 


2.643 


2.617 


96.090 


85.350 


89.292 


KBN 


2.685 


2.351 


2.620 


2.861 


85.652 


85.379 


70.497 


CaBN 


2.652 


2.304 


2.582 


2.839 


79.760 


87.374 


57.097 


GaBN 


2.558 


2.013 


2.617 


2.680 


77.240 


93.780 


43.431 


GeBN (Not K4 type) 


AsBN 


2.632 


1.922 


2.817 


2.684 


73.552 


99.772 


29.923 


SeBN 


2.629 


2.38 


2.610 


2.730 


93.992 


83.161 


89.371 


BrBN 


2.682 


2.361 


2.733 


2.739 


94.082 


85.632 


85.293 


RbBN (Not A'4 type) 


SrBN 


2.814 


2.446 


2.735 


3.015 


80.107 


87.151 


58.018 
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FIG. 10: (a) and (b) Distances and ratios of the nearest-neighbour distances between Ba and X (triangle), and X and Na 
(pentagon) for fully optimized BN, NaBN, MgBN, GaBN, FBN and CIBN crystal structures with K4 type frame within LDA. 
Classification of the B and N atoms (Ba and Na) is based on the irreducibility of the crystal symmetry as shown in Fig. 1. 



